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Abstract 

A construction of the heat kernel diagonal is considered as element of generalized Zeta 
function, that, being meromorfic function, its gradient at the origin defines determinant of 
a differential operator in a technique for regularizing quadratic path integral. Some classes 
of explicit expression in the case of finite-gap potential coefficient of the heat equation are 
constructed. 

1 Introduction 

A necessity of a Green function diagonal study is directly connected with the generalized zeta- 
function (GZF) theory of elliptic differential operators [1], which is successfully applied to a 
regularization of the operators determinants |2]. Such elliptic problems, for example, appear 
as Laplace transform of heat kernel equations, generally, with variable coefficients, which are 
conventially named as potentials. An important application of the theory is evaluation of 
semiclassical quantum corrections calculations to nontrivial classical solutions of important 
nonlinear equations and field theory [3l H]. The corrections are intimately linked to the fun- 
damental solutions of related linear problems for the heat operator Laplace transform, which 
diagonal enters the zeta-function definition. Such regularization, for example, is realized in 
explicit form for kink solutions of the integrable Sine-Gordon equation [3 Ej as well as non- 
integrable Landau- Ginzburg ((^^) models [H E]. The kink solution (as well as multikink one) 
in this context corresponds to the case of point spectrum of the elliptic operators that appear 
after division of variables. 

This paper is devoted to investigation of a wide class of potentials which spectrum is con- 
tinuous with eventual gaps - more precisely so-called finite-gap ones - see, e.g. the book [7]. 

Such potentials and, especially three-gap one, correspond to basic three-wave interaction, 
which is important in many quasiperiodic processes description, an exemplary applications one 
can find in [9]. 

In the Sec. 2, starting from Laplace transform of the heat equation by time, we derive 
a nonlinear equation for the Green function diagonal, along ideas similar ones, mentioned in 
[Bl E] in the context of other equations. We construct its solutions in cases which potentials 
has direct link to the polynomial functions in appropriate variables (Sec. 3). The last section 
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is devoted to examples and appendix contains a Mathematica program, related to a class of 
illustrations. 



2 The equation 

We are interested in a class of problems, connected with the parabolic partial differential oper- 
ator Green function (kernel of heat equation) 

d d"^ \ 

d^^d^" ^^^^ ) 9i^^^o,y) = S{x - xo)S{y), (2.1) 

where g{x, xq, y) & S is the fundamental solution over Schwartz space S; S{x — xo), S{y) are 
Dirac delta-functions. After Laplace transform: 

P+ - U{x) j g{x,Xo,p) = 6{x - Xo) (2.2) 

The construction of GZF in fact rely upon the Green function diagonal. In [B] a statement 
about g{p,x,x) = G{p,x), is used. Namely, G{p,x) solves the equation 

2GGI - {G'J' - 4{U{x) - p)G' + 1 = 0. (2.3) 

on condition, that U{x) is bounded. The equation resembles one of derived by Hermit for [8]. 
Proof: 

Let us consider homogeneous equation 

p+^-Uix)^ fip,x,Xo) = 0. (2.4) 

The fundamental solution of fl2.4p is built by standard procedure [?]. It has two linearly 
independent solutions, for example (j) and ip, converging respectively at — oo and +oo. One 
can represent through (p and i/j respectively for x < xq and x > xq with a sewing condition 
determined by equation (12. 2p 

gD{p,x,xo) = < {,) \ ^ . (2.5) 

^ ^ B[XQ)tp[p,x), X > Xq ^ ' 

From continuity condition oi (Jd one gets 

v4(xo)0(p,a;o) = B{xo)iJj{p,xo). 

What leads to: 

A{xo) = C{xo)i'{p,xo), 
B{xo) = C(xo)0(p,2;o). 

Due to the symmetry of Green function in respect to exchanging x and xq, C{xo) is constant 
(later referred as C). To obtain condition for derivatives of (p and ip one integrates (12.21) over x 
in an e neighbourhood of Xq: 

p+ - U{x)j go {p,x,xo) dx = 1, (2.6) 



dgp 
dx 



xo+e rXQ+e 

+ {p -U{x)) gD{p,x,Xo)dx = 1, 

x=xo—e Jxo-s 



dcj) dip 

— {p,XQ + e)Cilj{p,xo)-—{p,XQ-£)C4){p,xo)+ / {p-U{x))gD{p,x,XQ)dx = 1. 
dx dx J^^_^ 

In e — )■ limit above equation reduces to 

d(p dtp 

— {p, xo) Cip{p, xq) - -Q^{P, xo) C(j){p, xo) = 1. (2.7) 

Since solutions of fl2.4p are linear, one can assume C = 1. Then fl2.7l) reduces to 

d(p dip 

— (p,xo) ip[p,Xo) = —{p,Xo) p{p,Xo) + 1. (2.8) 

Actual proof will be made, by inserting (12. 5p to (12. 3p . For brevity function arguments will 
be omitted and ' will denote a derivative with respect to x 

2iPP {iP"(P + 2iP'p' + iP(P") - [ip'p + iPP'f - 4([/(x) - p)ip'^p'^ + 1 = 

2iP'^P {P" - {U{x) - p)(P) + 2iP4>'^ {iP" - {U{x) - p)ip) + AiP'p'iPP - [ip'p + iP(P'f + 1 = 0. 
Because of (12. 4p two first elements are nullified. One also uses property (12. 8^ : 

AiP'p'iPP- {2iP'(P+lf + 1 = 0, 

Aip'cp'ipp) - Aip''^(p^ - Aip'tp -1 + 1 = 0, 

iP'(P'iP(P - Vy^0^ - ip'(p = 0, (2.9) 
ip'^P"^ + ip'P - ip'^(p^ -^'0 = 0. (2. 10) 

Thus the proof is concluded. 

It is important to note, that the transition is general and doesn't rely on the nature of U{x) 
as long as it's bound. It's usefulness is dependent on a few qualities of the potential though. 

3 The main equation solution 
3.1 Substitutions 

We consider a class of solutions of the equation (12. 3p . to be written in a form 

It is most useful, if there exists a variable transition x ^ z, U{x) ^ u{z), in which P and Q 
are polynomials. Basic conditions for it to be possible are: 

1. n is a polynomial in z, 

2. {z'^y is a polynomial in z (note, that z'J.^ = ^^(-z^.)^)) 



This does not ensure simplicity of solutions as will be shown further in the text. At this point, it 
is important to notice, that the second condition restricts z{x) - apart from a class of elementary 
functions - to elliptic and hyperelliptic functions - see, e.g. [7]. 

Let us assume, that {z'^)"^ and u are polynomials in the variable z of degree L + 1, and K 
respectively, hence, polynomial of the degree L. We also assume, that coefficient by 

the highest power term of {z'^)^ is equal to 1 (which is always attainable). After the change of 
variables, the equation will take the form: 

2P(P"(4)2 + - {P'z'J - 4(n(.) - p)P' + 4Q = (3.12) 

We also assume solution in a given form: 



QiP) = E'=0 ' QnP- 



(3.13) 



3.2 Classification 



We will now proceed to analyse the solution by separating the equation in respect to powers of 
p and z. The equation for p^, ^2Mo+max(x,L-i) ^g^j^gg following form: 
If > L - 1 

4mxPo,Mo = (3.14) 

UK <L-1 A Mo > 1 

2Po'mo(Mo(Mo - 1) + ^Mo) - Po%/„Mo^ - ^ukPImo^k^l^i = (3-15) 

+ {L- l)Mo - 4ukSk,l~i = (3.16) 

Note, that Mq = leads to K = (this case will be examined later in the text). Another 
conclusion is, that < L — 1 is necessary for sought type of solutions. Furthermore, this leads 
to following, more precise conditions: 

L > 1 {due to K>0) . . 

K = L — 1 if or Mq to have positive value) 



Equation for p^A^+i 



4(5^^^.'^') +4g2iVH-i = (3.18) 



leads to following conclusions: 

Mn = a Plo = -q2N+i- (3.19) 
Let us now look at subsequent equations for descending powers of p. For p"^^ we have 



Mjv-i 

Mz)Pn,o + 8Pn,o Y1 + 4g27v = 0, (3.20) 

/=0 



which leads to 

Mn-1 ^ ^ 

E ^^-1.'^' = 2 - l^j ' (3.21) 

M^„i = (3.22) 
For p2Af-i gg|^ (^fQj, ^]^g highest power of 2): on condition Mn-2 > ^Af-i + K {z ^ x"^) 

4:Pn,oPn-2,Mm^2 — 0' (3.23) 

on condition Miy^2 < Mjy^i + K 

2Pn,oPn-i,m^.AMn-i{Mn-i - 1) + ^M^_i) 

— 2F/v,o-P/v-i,A/jv_i — 8ukPn,oPn-i,Mn-i (3.24) 

+8P/V,0-P/V-2,A/Ar_2^A/]V-2,-J--l + 4q'2Ar-l = 0. 

It's obvious, that M7V-2 < ^tv-i + -ftT is a necessary condition for (jlj). Now we can consider a 
general rule for all remaining equations. Thesis: \/o<k<N-iMk < {N — k)K. Proof by induction: 
If Mk > {N — k)K, then equation for p^+'^+i and highest power of z takes form: 

4P7v,o^fc,A4 = (3.25) 

If < {N — k)K and Wk<i<NMi = {N — l)K (possibility giving the highest possible value of 
Mfc), then equation for p^+'^+i and highest power of z takes form: 

2P^,oPfe+i(Mfc+i(Mfc+i - 1) + ^Mfc+i) 

+'^Yl,n=k+2Pri,MuPN-n+k+l,MN-n+k+ii^N~n+k+l{MN^n+k+l " 1) 
-n+fc+l) ~ 2 'Ylin=k+2 -Pn,M„ -P/V-n+fc+l^TV-n+fc+l 
-4'Ui<: Z]n=fc+1 Pn,MnPN-n+k+l 



(3.26) 



+4 Z]n=fc+1 Pn,MnPN-n+k,MN~rv+k 
,Mk+i+L—lPNfiPk,Mi^ 



Thus the solution exists only if the thesis holds. This leads directly to a minimal condition on 
N: 

N>^ VMo > K (3.27) 
K 

In summary: existence of solutions of form (jlj) depends on the power of the potential (K), 
power of {z'^)"^ (L < 1 can give abnormal results), amplitude of the highest power term of the 
potential and there exists a definite formula for the minimal value of N. 



3.3 Solving algorithm 

1. If there exists a solution in the form (jl]) for a given (which fulfils requirement f l3.3p ). 
it can be obtained in a straightforward manner. Since the actual value of is un- 
known, one starts with the minimal possible value After separating the equation 
in respect to powers of p one analyzes the resulting equations starting from the highest 
power of p. All those equations can be written in a manner similar to fl3.2ip (here for the 
p'^^'^^^withpossiblenfromN-ltoO): 

{N~n)*K 

Pn,o Y1 Pn,iz' = F{PN,o,PN-i,K,PN-i,K-i,---,Pn + l,0,z)-^^^^ (3.28) 

1=0 



where F contains all elements of equation not written explicitly. It is easy to see, we can 
obtain all coefficients, except for Pn,o, as solutions of linear equations, since all elements on 
the RHS, except for qN+n+i are known. As for the equation for z^, it is more convenient, 
to express qN+n+i in terms of P„,o- Solving all equations down to p^+i gives us all Pn,i as 
well as some of g„ in terms of {Pn.o}ie{o,...,N}- It is important to note, that all calculations 
done up to that point stay relevant, even if value of will have to be increased. 

2. In the next step, we use the equation for to calculate the possible values of Pn,o- Again, 
if we start from the highest powers of z, we can obtain those coefficients as solutions of 
linear equations, since any element containing P„ q is proportional to at most z^^^~'^~^^^ 
and any element containing P^q is proportional to at most (negative exponent 
means, that such coefficients are not present in equation for p^). 

3. Subsequently we check the solution. If it doesn't hold, we increase the value of n, add 
relevant components to P and Q polynomials, calculate values of all new coefficients and 
go back to step 2. 

Since g„ aren't necessary for calculation of any Pn^i and checking the solution, we can slightly 
simplify the algorithm by calculating g„ after all other coefficients. 

Described algorithm was implemented in Mathematica 7 and used to obtain solutions presented 
in section m 

3.4 On uniqueness of solutions 

Using the above algorithm we obtain all coefficients as solutions of linear equations in respect 
to sought coefficients. This leads to a simple conclusion, that for a given A^ solutions of form 
(jl]) are unique except for constant u, in which case the solution has no dependence on z and 
because of this, there is no relation between any of Pnfi (14.11) . As yet, there is no method of 
finding all allowed A^ for a given potential. Therefore uniqueness of solutions is uncertain. 



4 Exemplary solutions 

4.1 Constant potential 

Let's consider a constant potential 

U{x) = u (4.29) 

It is obvious, that no change of variables is necessary, thus we can use z = x {L = —1). 
Equation for p^z"^^^" immediately gives 

- iuPl,,^ + 4qo6Mo,o = (4.30) 

Since u can have an arbitrary value, if we shift the p variable, this equation only holds for 
Mq = 0. This means, that the simplest solution would be 



(4.31) 



This is not the only one, as will be shown. Let us consider a solution for potential fl4.29p and 
an arbitrary N. Equation for p^A^+i gjyes as usual 

Mtv = (4.32) 

q2N+i = -Pi, (4.33) 



Mn^i = (4.34) 
Pm-1,0 = l (uPn,o - ^) (4.35) 



Equation for p gives 



2 V ' Pn,o, 

Subsequent equations will look alike, with Mj = for all i. It is easy to see, that one obtains a 
total of 3A^ + 3 parameters with only 2N + 2 equations and one can obtain a solutions for any 
value of A^. In a sense, it's a consequence of the condition fl3.27p . since ^ is an indeterminate 
symbol. 

4.2 Triple-gap cnoidal potential 

Let's take a solution one order higher then that for 0^ cnoidal solution: 

U{x) = -Um^k'^ cn^{mx; k) (4.36) 

z = cn^{mx] k) (4.37) 

{z'^f = Am^z{l- z){l-k^ + k^z) (4.38) 

4'^ = 2m2(-3A;2^2 + (4A;2 -2)z + l- k"^) (4.39) 

Mo = 3 (4.40) 

iV = 3 (4.41) 
Algorithm explained in section (13. 3p gives (assuming P30 = 1 for simplicity): 

P2 = -2m2(7 + /fc2(-14 + 3z)) (4.42) 

Pi = 771^^(49 + A;2(-256 + 78z) + A;'^(256 + 3z(-52 + 15z))) (4.43) 

Pq = _3^„6(^i2 + 8A;2(-19 + 9z) + 3^^(128 + z(-121 + 45z))+ (4.44) 

fc^(-256 + 3z(121 + 5z(-18 + 5z)))) (4.45) 

Q{p) = -([-4 + 8k'^)m'^ +p){{9-96k^ + 96k^)m^ + 10{-l + 2k'^)m'^p + p^) (4.46) 

((9 - 42A;2 + 33A:^)m^ + 2(-5 + 7k^)m^p + /) (4.47) 

(3A;2(-8 + llk^)m'^ + 2(-2 + 7k^)rr?p + p^) (4.48) 



5 Conclusion 

Described equation allows calculation of heat equation's Green function diagonal in a straight- 
forward manner. Developed algorithm should be especially useful for finding solutions for 
finite-gap potentials, which naturally emerge in periodic and quasi-periodic structures. 
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